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ABSTRACT 

We report the results of a spectroscopic study of the high-mass protostellar object NGC 7538 IRS 9 
and compare our observations to pubHshed data on the nearby object NGC 7538 IRS 1. Both objects 
originated in the same molecular cloud and appear to be at different points in their evolutionary histo- 
ries, offering an unusual opportunity to study the temporal evolution of envelope chemistry in objects 
sharing a presumably identical starting composition. Observations were made with the Texas Echelon 
Cross Echelle Spectrograph (TEXES), a sensitive, high spectral resolution (i? = A/AA ~ 100,000) 
mid-infrared grating spectrometer. Forty-six individual lines in vibrational modes of the molecules 
C2H2, CH4, HCN, NH3 and CO were detected, including two isotopologues (^^CO, ^'^C^^O) and one 
combination mode (1/4 -|- C2H2). Fitting synthetic spectra to the data yielded the Doppler shift, 
excitation temperature, Doppler h parameter, column density and covering factor for each molecule 
observed; we also computed column density upper limits for lines and species not detected, such as 
HNCO and OCS. We find differences among spectra of the two objects likely attributable to their 
differing radiation and thermal environments. Temperatures and column densities for the two objects 
are generally consistent, while the larger line widths toward IRS 9 result in less saturated lines than 
those toward IRS 1. Finally, we compute an upper limit on the size of the continuum-emitting region 
(^2000 AU) and use this constraint and our spectroscopy results to construct a schematic model of 
IRS 9. 

Subject headings: stars: formation — stars: pre-main sequence — techniques: spectroscopic — ISM: 
individual objects (NGC 7538 IRS 9, NGC 7538 IRS 1) 



1. INTRODUCTION 



Details of the pre-Main Sequence (pre-MS) life cycles 
of low-mass stars such as the Sun [M ^ IMq) have been 



gleane d from decades of research; see |McKee fc Ostriker] 
( 2007 1 for a recent review. This is thanks in part to 
the relative amenability of low-mass protostellar objects 
to study. However, details of high-mass {M > 8Mq) 
star formation are not as well understood at present, in 
part because high-mass stars are not generally forming 
near the solar neighborhood and remain deeply embed- 
ded during their pre-MS evolution. An important unre- 
solved issue is whether high-mass stars form in a man- 
ner rese mbling a scaled- up ve r sion of low-mass star for- 



mation (Keto & Zhang 2010 Johnston et al 



2011 



through other mechanisms uncharact eristic of the 



or 
ow- 



mass case such as coll isions /m ergers fBonnell et al.|1998 
Bonnell fc Bate|2005[ [Bally fc Zinnccker 2005' ) and com 



petitive accretion m clusters ( |Bonnell et al.,|2001 2004 ) . 
Clearly, the environments in which high-niass star tor- 
mation takes place seem to influence the formation pro- 
cess itself. Massive stars are most often seen to form in 
large groups ( "OB associations" ) near the e dges of dense 
clumps in g ia nt molecular clouds (GMC s; Zinnecker & 



Yorke|[2007|. [Elmegreen fc Lada| ( |1977[ ) first suggested 
that spatially distinct subgroups of stars within OB asso- 
ciations trigger waves of self-propagating star formation 
in GMCs via shock/ionization fronts that progress from 
the outside in. This process may account for observations 
of OB subgroups that lie in linear arrangements along 
the Galactic plan e, showing a monot onically-increasing 
sequence of ages (Blaauw 1964 1991). Rapid improve- 



ments in both observational capabilities and theoretical 
sophistication in recent years are advancing our knowl- 
edge of the circumstances of massive star formation and 
enabling these fundamental questions to be addressed. 

The origin of these stars is of particular interest due 
to their role in the Galactic ecosystem. In addition to 
serving as the principal sites of heavy element nucleosyn- 
thesis, the violent end of massive stars in supernovae in- 
jects considerable mechanical energy into the inte rstellar 
medium (ISM) while enriching it with metals (Arnett 
| 1996[ ). The debris are reconstituted in molecular clouds 
from which new stars form, often accompanied by cir- 
cumstellar disks ( Jiang et al.||2008 ). These disks are ev- 
idently the source of accretion onto high-mass stars as 
they approach the MS in the Hertzsprung-Russell dia- 
gram and ignite hydrogen burning, then continue t o ac- 
crete their f inal masses, moving parallel to the MS ( Zin- 



necker fc Yorke__2007J . Relatively little is known about 
the structure and chemical composition of massive proto- 
stellar envelopes, and even less about how they evolve in 
time. Ready comparison with envelope models for low- 
mass objects is instructive but only of limited use, in part 
due to the radically different radiation environment in 
the high-mass case. The composition and structure of the 
envelope are of interest as they determine the final mass 
of the star through an accretion mechanism whose details 
are still unclear. Piecing together the components of pro- 
tostellar envelopes is complicated by the observation that 
most objects are deeply embedded, particularly early in 
the formation process. Spectroscopy of molecules is use- 
ful in delineating physical structures at high extinction 
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where other methods are Umited. Pure rotational tran- 
sitions of molecules in the millimeter probe a sufficiently 
low energy regime but such observations are typically of 
low spatial resolution. Better resolution is achieved in 
the infrared where many ro-vibrational molecular lines 
are located. 

We present high resolution, mid-infrared spectra of a 
number of molecules toward an embedded infrared source 
in the NGC 7538 starforming complex, NGC 7538 IRS 9, 
and compare the results to a previous s tudy of the nearby 
source NGC 7538 IRS 1 (| Knez et al.|2009) . Both sources 
are presumed to harbor high-mass protostellar objects in 
their interiors. Their proximity to one another suggests 
that, having formed within the same molecular cloud, 
the starting chemical abundances for each object should 
have been substantially similar. A comparison of their 
current states provides insights about differences in their 
ages and/or divergence in their individual evolutionary 
histories. 



Weak fre e-free emission (< 60 /iJy at 3.6 cm; Sandell 
et al.|[2005 l is observed toward IRS 9. Further evidence 



tor the comparatively primitive nature of the envelope of 
IRS 9 is given by the detection of various ices in its In 



frared Space Observato ry (ISO) spectrum (Whittet et al. 

These ices are not seen toward 



Gibb et al. 2004) 



1996 

IRS 1 j which appears to be in a relatively advanced point 
in its pre-MS evolution. It is a considerably stronger 
source at centimeter wavelengths and shows evidence of 
both ionized and molecular components in a bipol ar out 
flow originating in a disk seen nearly edge-on (Sandell 
et al.||2009||Qiu et al.l|2011[ ). IRS 1 also shows signs of a 
circumstelJar disk and an emerging compact H II region; 
it is thought to be a forming late O star undergoing active 
accretion, s uggestin g a more advanced ev ol utionary state 



( Lacy fc Jaffc 2007 
2009| further 1^ 



Beuther et al. 



Sandell et al. 



2012[ ) 

hat the system is heavily accreting 
(M ~ 2 X 10"'* Mq yr~^) and that the accretion may 
be episo dic in nature . On the basis of GHz line obser- 
vations, Surcis et al. (20111 conclude that the accretion 
appears to be the result or radial infall from a "torus" 
thousands of Astronomical Units (AU) in size, traced by 
H2O and CH3OH masers, instead of the Keplerian disk 
favored by Pestalozzi et al. ( 2004 ) . 

Until recently, estima tes of the distance t o NGC 7538 
placed it at -3.0 kpc (iHeyer et al.||2001j [Brunt ||2003 
Foster & Routlcdgc 2003; Balog et al. 2004VRcidirWir 



son^2005^,Kameya. 2006^ .Araya^et al.^.20(37J. Ho wever, we 



adopted the distance oflMoscadelli et al. } 2009 ) , who find 



a value of 2.65^q kpc from trigonometric parallax mea- 
suremen ts. The total luminosity of IRS 9 is about 3.5 x 
10'' Lq ( [Sandeh et al.||2005t corrected to the Moscadelh 
et al. dist ance) and it is a brigh t IRAS 12 /xm source 
at ~60 Jy ( Beichman et al. 1988 ) . A zero-age Main Se- 



quence (ZAiVlSj star with 



ais luminosity would have a 



spectral t ype of B O. 5, an effective temperature T^t f = 
26200 K | Panagial973 |, and a mass o f ~13 Mq (|Et 



str5m et al. 2012 ). Boogert et al. (2004) suggest the ex- 
istence of an inner, warm molecular emission region with 
a radius of ~70 AU; the dust temperature a tthis dis- 
tance for an object of IRS 9's luminosity would be ~1000 
K depending on the dust mass opacity. Millimeter data 
show that IRS 9 is a site of active, massive star forma- 
tion ( Sandell et al.|[2005 ) in which young stellar objects 



(YSOs) are driving a set of bipolar molecular outflows on 
a dynamical timescale of < 20,000 years. At lea st one of 
these outflows app ears to be highly energetic. [Mitchell 
& Hasegawa ( 1991 ) report the discovery of an outflow in 
IRS 9 with a velocity of 110 km s~^, which they interpret 
to have emerged as recently as 1,200 years ago, perhaps 
indicating the beginning of a recent episode of accretion. 
Observations of these prot ostellar ob j ects a re broadly 



consistent with the Elmcgrcen & Lada (1977) picture of 
sequential high-mass star formation, as elabor ated on 
by others hi the particular cas e of NGC 7538 dWernerl 



et al.[[1979| [Fischer et a l 



Campbell & I'hompson 



1980 Dickel et al. 1981). The 



ggymodel of NGC 7538 holds 



that a hrst generation of massive stars caused the visible 
H II region, which expanded and impinged on the neigh- 
boring molecular cloud; this compressed the cloud and 
generated shocks that caused the observed near-infrared 
fluorescent H2 emission northwest of IRS 1. In this pic- 
ture, high-mass star formation is progressing toward the 
southeast, and IRS 9 belongs to a second wave of star 
formation as it is placed further along the direction of 
propagation of the shock front. These objects then rep- 
resent two snapshots of the evolution of a massive young 
protostellar object in a particular cloud and are useful in 
addressing temporal evolution questions. 

This paper is organized as follows. In Section 2 we 
give details concerning how we obtained and reduced the 
data, followed by a detailed description in Section 3 of 
the molecular species detected. The data analysis is de- 
scribed in Section 4 with emphasis on spectral modeling, 
and we give some interpretation to the results in Section 
5. We summarize these results and draw some conclu- 
sions in Section 6. 

2. OBSERVATIONS AND REDUCTIONS 
We observed NGC 7538 IRS 9 with the Texas Eche- 



lon C ross (X) Echelle Spectrograph (TEXES; Lacy et al. 
2002[ ) on the Gemini North 8 m telescope in 2007 Oc- 
tober. Similar observations of IRS 1 were carried out 
with TEXES at the NASA Infrared Telescope Facihty 
(IRTF) 3 m telescope in 2001 June, 2001 November, 2002 
September, 2002 December and 2005 November. Cir- 
cumstances for all IRS 9 observations are given in Ta- 
ble [H I RTF observations o f IRS 1 form the basis for the 
work of Knez et al. ( 2009 ) . Additional observations of 
IRS 1 were made at Gemini in 2007 October to supple- 
ment previous TEXES data taken at the IRTF. 

TEXES is a cross-dispersed grating spectrograph de- 
signed to operate at high spectral resolution (R = 75,000- 
100,000) over a range of mid-infrared wavelengths from 
5 to 25 /im. Light diffracts off a coarsely-ruled "eche- 
lon" grating (jMichclson 1898) at a large incidence angle 
with respect to normal at order n ^ 1500. The light 
then diffracts off a cross-dispersion grating and is im- 
aged onto a 256 x 256 pixel Si:As detector as a series 
of orders. At a given spectral setting between five and 
ten cross-dispersed orders are imaged onto the detector 
for a spectral coverage of about 0.5%. The continuity 
of orders is wavelength-dependent; below 11 /im order 
overlap is sufficient to yield continuous coverage over an 
entire spectral setting. At wavelengths longer than 11 
fim there are gaps between adjacent orders. 

On Gemini the pixel scale of TEXES was --0'.'14 while 
at the IRTF it was ~0"36. The pixel scales also compare 
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TABLE 1 

Circumstances of TEXES observations of NGC 7538 IRS 9 



Wavenumber 

UT Date Telescope Center (cm~^) 



2001 June 15 IRTF 768 

780 

2001 June 18 IRTF 744 

2001 June 26 IRTF 761 

2001 June 28-29 IRTF 734 
2007 October 21 Gemini 828 

1308 

2007 October 25 Gemini 930 

743 

2007 October 26 Gemini 853 

1323 

2007 October 27 Gemini 2055 

2007 October 28 Gemini 730 
2007 October 30 Gemini 2085 



Included Features 



C2H2 -R(15,16) 
[Nell] (2Pi/2 ^ 2P3/2), C2H2 U5 

i?(21,22), HCN 1/2 iJ(22,23) 
C2H2 V5 -R(5,6), HCN U2 i?(10) 
C2H2 !/5 ^?(12,13), HCN U2 R{15,16), 
HD, H^^CN U2 -R(18) 
C2H2 1^5 HCN V2 H(6,7) 

HNCO VA P-branch, ^-^NHs V2 sP(7,K) 
12CH4 Vi R{0), 13CH4 U4 R{1) 
"NHg 1^2 aQ(J,K) 

C2H2 U5 P{5,6) 
i^NHa 1/2 aP(4,K) 
12CH4 U4 R{2), 13CH4 U4 R(4), 
C2H2 U4, + us P{l-3) 
OCS U2 P(12-22), i^co v=l-0 P(21,22), 
"CO P(ll,12), Ci«0 P(10,ll), C"0 P(16) 
C2H2 Q-branch 
i2cO v=l-0 P{14-15), "CO P{3,4), 
C"0 P(2,3), Ci'^O P{8,9) 



favorably with the diffraction hmit at 10 /im of each tele- 
scope, ff.'TQ and 0'.'26 at the IRTF and Gemini, respec- 
tively. In both cases the seeing disc is, at worst, critically 
sampled under good seeing conditions. All observations 
were made at a resolving power R ~ 80,000, or Av ~ 
3-4 km s~^ and the spectral coverage at each spectral 
setting was AA « 0.06 /im. The spectral sampling of the 
TEXES array is about 1.0 km s'^ Slit widths were O'.'S 
at Gemini and 1'.'4 at the IRTF. The slit length was ~4" 
on Gemini and ~10" on the IRTF, varying with wave- 
length. We alternated the telescope pointing between 
two positions on the sky in the slit direction, separated 
by ^ 1-5 , during an exposure sequence. Subtraction 
of adjacent "nod" positions removes the contribution to 
the signal from night sky emission. The separation of 
points in the nod, or "nod throw" , is determined by the 
slit length. 

To correct for telluric atmospheric absorption we ob- 
served comparison objects whose spectra are close to 
blackbody over the wavelength range of observation. At 
wavelengths beyond ^ 8 /im, asteroids meet this require- 
ment while providing the most flux of available object 
types. Toward 5 /tm, most asteroids provide insufhcient 
flux to achieve a high signal-to-noise (S/N) ratio, so we 
used bright, hot Main Sequence stars such as Vega (a 
Lyr), Fomalhaut (a PsA) and Mirfak (a Per) as divi- 
sors. We aimed to observe the divisor for a given target 
at nearly the same airmass as the target observation for 
better divisions. Differences in airmass at the times of 
observations may be corrected after data reduction by 
taking the divisor spectrum to the power of the ratio of 
the airmasses. Residual structure in divided spectra was 
removed by fitting and subtracting low-order polynomi- 
als. 

Finally, given the high spatial resolution of the Gemini 
observations, we considered the possibility that TEXES 
might resolve the emitting/ absorbing region around IRS 
9 if it were spatially extended. We obtained a series of 
spatial-spectral scan maps of IRS 9 that allow us to place 
an upper limit on the linear size of any resolved struc- 



ture. The approach and results are discussed in detail in 
Section 15.51 

2.1. Data Reduction 

Our data reduc tion pipeline is described in detail in 
Lacy et al. ( 2002[ ). Flux calibration is based on the ap- 
proach of calibrat ing millimeter and su bmillimeter data 
in the manner of Ulich & Haas (19761. A typical sky 
calibration cycle involves collecting data on the sky and 
object with two positions of a rotating chopper blade 
mounted just above the entrance window to the TEXES 
Dewar: one painted flat black and one low-emissivity, 
polished aluminum ("shiny") surface. The result of di- 
viding by the black-sky difference removes atmospheric 
absorption and serves as both a spatial and spectral flat 
fleld image provided that the black position, sky, and 
telescope temperature s are approximately equal. We 
quote the Lacy et al. (2002) expression for the object 
intensity as a function ot frequency, 



Iu{ohject) « Sv{object — sky) 



S^{black — sky) 



(1) 



where S^, is the measured signal and is the black- 
body function at the telescope temperature Ttei- In the 
approximation that /j, (black)-/^ (sky) — B^{Ttei)e~'^''*'y ^ 
this procedure corrects for sky absorption. Wavelength 
calibration is achieved by using the known wavelengths 
of night sky emission lines; expressed in velocity units, 
this method is accurate to ^1 km s~^. 

Pipeline software written in FORTRAN performs the 
reduction steps resulting in a wavelength- and flux- 
calibrated object spectrum along with an estimate of the 
noise, the sky transmission spectrum, and the 2-D unex- 
tracted slit image. The program also corrects for optical 
distortion, removes cosmic ray spikes and fixes bad pix- 
els. The systematic accuracy of radiometric calibrations 
is ^-^25%, ignoring telescope pointing uncertainties. 

3. DESCRIPTION OF THE DATA 
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We observed IRS 9 in the lines of seven molecules at 
the settings summarized in Table [l] The selection of set- 
tings was motivated by the desire to include as many 
species and individual lines as possible given the severity 
of telluric atmospheric absorption at a particular setting. 
In particular, variations in the II2O column over Mauna 
Kea often influenced the settings observed on a partic- 
ular night, as many of the settings are sensitive to the 
intensities of telluric water lines. The systemic velocity 
of IRS 9 and motion of the Earth with respect to the 
Local Standard of Rest (LSR) were taken into account 
in choosing settings in which certain lines were shifted 
either into or away from nearby telluric absorption fea- 
tures, and in many cases required trade-offs in order to 
maximize efficiency and scientific return. Details of the 
observations for specific molecules follow. 

3.1. C2H2 

At Gemini, two lines in the acetylene 1^5 symmetric 
bending mode R branch were observed (J— 5,6) along 
with approximately ten lines in the Q branch. The R- 
branch lines are shown in Figure [Tl while Figure [2] shows 
the Q-branch spectrum along with a simple LTEnt at a 
single rotational temperature of 100 K. This is discussed 
further in Section 4. Each acetylene line appears to con- 
sist of a single component. The -R(5) line is stronger 
than R{6), as is expected due to its three-times-greater 
nuclear statistical weight. i?(6) appears to have a weak 
redshifted emission component at an LSR velocity of -54 
km s~^ although this effect may be due to a poor fit to 
the continuum there. This component can be more easily 
seen in the top-center panel of Figure [lO] Earlier IRTF 
observations of the 1^5 R{^) hne (wavenumber center 734 
cm^^) in 2001 June revealed a positive detection, but 
more importantly, higher J settings i?(15,16) (768 cm^^) 
and i?(21,22) (780 cm^^) were non-detections. The lat- 
ter setting included the ^^3/2 fine structure line 
of [Ne II] ; consistent with the weak free-free emission ob- 
served toward IRS 9, this line was als o a non-detection 
and is discussed in detail in Section|44j Finally, although 
the positions of lines of the i3q12(jjj^ isotopologue fell 
in some of our spectral settings, we did not detect any. 

We also did not detect any lines in the 2^5 — ly^ "hot 
bands" of acetylene toward IRS 9. However, the hot- 
band lines in the spectral settings observed were not par- 
ticularly favorable for detection, so we cannot draw any 
strong conclusions about the 1/5 population from these 
observations. The observation of 1^4 -I- combination 
mode lines, shown in Figure [3] is of particular inter- 
est. This is because it allows a measurement of the 
C2H2 abundance with intrinsically weaker lines at dif- 
ferent wavelengths from the z^s band lines. Since the 
+ hne strengths are ~5 times weaker than the 1/^ 
line strengths, absorption in the 1^4 + 1^5 lines may be 
dominated by high column density clumps of gas. 

An energy level diagram indicating relevant vibrational 
levels of C2H2 is shown in Figure Hj The mode 
is the symmetric bending mode in which the H atoms 
move together. In the V4 antisymmetric bending mode, 
they move in opposition, resulting in no dipole moment. 
Therefore, transitions that change the 1/5 quantum num- 
ber are relatively strong, with Aui ~ 10 s~^. These 
include the v^, 21/5 — (between 2i^5 and 1^5 ), and 
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Fig. 1.— TEXES spectrum toward NGC 7538 IRS 9 showing 
features due to C2H2 and HON. From left to right, the asterisks 
mark the positions of C2H2 1^5 R{5), HON U2 -R(IO), and C2H2 
Us R{6). The upper trace shows the data, while the lower curve 
indicates the relative atmospheric transmission on the same scale. 
The HCN U2 -R(IO) line is a detection despite the telluric lines at 
that wavenumber not completely dividing out. Gaps in the spec- 
trum occur between echelon grating orders and where the telluric 
absorption is too great to be divided out. 

+ 1^5 — V4 (between 1/4 -\- vc, and V4) transitions. Tran- 
sitions that only change V4 do not occur radiatively, but 
transitions changing both V4 and 1^5, including V4+u^ and 

— V4, do occur although they are ~ 10 times weaker 
than those which only change v^. Because of these se- 
lection rules, 7.6 ^m photons can be absorbed in 1/4 -\- 
lines, but molecules excited into the V4 + 1/5 level most of- 
ten decay to the V4 level from which they can only decay 
to the ground vibrational level collisionally. 

The ratio of the V4 -I- P{i) and P{2) lines is difhcult 
to determine due to poor telluric division, but appears to 
be less than the 9:2 ratio expected from their rotational 
and nuclear statistical weights. This probably indicates 
saturation of the lines in spite of their intrinsic weakness 
relative to the transitions. We also note that the spec- 
tral region containing the 1^4 + lines is contaminated by 
dense forests of HNO3 and SO2 lines, giving the super- 
ficial appearance of bad flat-fielding. Another difference 
between the and the V4 + transitions involves the 
shapes of the lines; the V4 -\- lines are more asymmet- 
ric, with blue wings extending to nearly -90 km In 
this regard they resemble the P Cygni-like shapes of th e 
CO lines we observed toward IRS 9 (see Section 3.5). 
The 1^4 -f lines may result from absorption in a high 
column density component with a small covering factor 
that is strongly saturated. Or they may be formed in 
highly dust-obscured gas that is more visible at 7.6 than 
13 /im due to the lower dust opacity off of the silicate ab- 
sorption feature. The presence of the similar -80 km s~^ 
absorption feature in the CO lines supports this possi- 
bility. Alternately, the V4 -f vc, line strengths relative to 
the I's lines may be affected by oth er radiative transfer 
effects, as we discuss in Section 15.3 



3.2. HCN 

Three of our spectral settings contain i?-branch lines 
of the V2 bending mode of HCN with centers at 734 
cm-i (6 < J < 7), 745 cm"! (J=10), and 761 cm"! 
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Fig. 2. — Spectrum of the 1/5 Q branch of C2H2 observed with TEXES toward NGC 7538 IRS 9 (thin hnes). A fit is superimposed from 



an LTE pure absorption model with T = 100 K, 6 = 1 km s~^, and 7V(C2H2) 
been convolved with the TEXES instrumental line shape function. 



2 X 10 cm (thick line). The synthetic spectrum has 



+ 




0.95 - 



0.90 



Fig. 3. — Spectra of two P-branch lines of the 1^4 + combi- 
nation mode of C2H2 shown on an LSR velocity scale. The two 
spectra have been offset vertically for clarity, and the continuua 
are indicated with heavy straight lines. Both spectra are affected 
by strong, poorly-corrected telluric absorption at small blueshifts. 

(J=16). The i?(6) line fell into a region where atmo- 
spheric transmission is not favorable, and R{1) appears 
to be a marginal detection because of a weaker, nearby 
atmospheric feature. i?(10) (Figure IlJ is a probable de- 
tection. i?(16) falls close to the i?(13jline of C2H2 but is 
sufficiently separated from the acetylene line to prevent 
blending and is the only clean HON line detected. The 
i?(15) line of HON fell just outside the first order of the 



761 cm~^ setting. The S/N of the HON observations is 
generally lower than that of C2H2 R branch, but the 
HCN lines are clearly the broader of the two. The mean 
full width at half maximum (FWHM) of the four C2H2 
i?-branch lines is 6.5±1 km s^^, whereas the HCN i?(16) 
fine FWHM is ~15 km s^^. Finally, lines of the Hi^CN 
isotopologue are also found in these settings but not de- 
tected: i?(9) in 734 cm"! and i?(18) in 761 cm"!. i?(18) 
falls into a gap between orders and i?(9) was affected by 
poor atmospheric transmission. 

3.3. CHi 

We obtained data on three i?-branch lines of CH4: 
the i?(0) line at 1311.43 cm"-'^ and the two i?(2) lines at 
1322.083 cm-i and 1322.152 cm"!. In addition, the set- 
tings for these lines contained the positions of the i?(l) 
and i?(4) lines of ^'^CH4, but both of these lines were clear 
non-detections. The ^^CH4 lines show clear asymmetries 
indicative of multiple velocity components. Each feature 
has a principal component corresponding to greatest ab- 
sorption at Vlsr — -65.0 km s^^ with additional compo- 
nents on the red wing of this feature separated in velocity 
by up to 6.5 km s^^. The spectra are plotted in Figure[5| 

3.4. NE-i 

Two spectral settings observed contain a total of ten 
P-branch lines of ammonia, three in the V2 "umbrella" 
mode antisymmetric P(4,X) group and seven in the vi 
symmetric P{7,K) group. Of the possible lines in these 
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Z/4 + i/5 



13.7 fim 



13.7 fim 



Fig. 4. — Energy level diagram for some mid-infrared-active, ro- 
vibrational transitions of C2H2. Radiatively-allowcd transitions 
are denoted with solid lines and coUisionally-allowed-only transi- 
tions witli dashed lines. The excited vibrational levels are split 
by interactions between the rotational angular momentum and the 
angular momentum of the bending-mode vibration; however, the 
energy differences between these levels are too small to be seen at 
the resolution of the diagram. 




-80 -60 

Vlsr C^fTi s"') 

Fig. 5. — Spectra of three features in the 1/4 band of CH4 toward 
IRS 9, shown on an LSR velocity scale. For the spectrum contain- 
ing the R{2) lines, the velocity scale has been set to the LSR value 
for the stronger component (rest frequency 1322.085 cm~^). The 
individual spectra have been offset vertically for clarity. 

settings, only aP(4,2) was not observed because of its 
proximity to a strong telluric atmospheric feature. Al- 
most all of the NH3 lines share a characteristic P Cygni- 
like profile, with an emission peak near the systemic ve- 
locity, around -60 km s~^, and weak blueshifted absorp- 
tion separated in velocity by 2.5 to 10 km s~^. This may 
be compared to the -56.5 and -59.8 km s~^ absorption 
components seen in mid-infrared NH3 absorption transi- 



tions toward IRS 1 ( jKnez et al.|2009 ) and a value of -60.0 
km s~^ for the (3,3) rotational transition of NH3 in Very 
Large Array (VLA) obse rvations of IRS 1 at centimeter 
wavelengths ( [Henkel et a l. 1984). The degree of absorp- 
tion toward IRS 9 relative to emission varies as well; the 
highest K lines in the sP{7,K) series show no absorp- 
tion at all despite the fair S/N of the data. Finally, we 
observed a setting near 930 cm~^ containing the posi- 
tions of the 1/2 aQ-branch lines (1 < J < 9); the data 
show weak absorption features at some of the indicated 
points in Figure [6j For comparison, a similar plot of the 
sP{7,K) branch spectrum is also shown in the figure. We 
did not observe lines in the 1^2 i?-branch of NII3 because 
of their proximity to regions of significant telluric atmo- 
spheric absorption. Also, the i?-branch lines are located 
in a spectral region of rapidly increasing absorption due 
to the 9.7 fim silicate dust feature. This fact significantly 
impacts our interpretation of the NII3 P- and i?-branch 
lines, as we discuss in Section [53] 



3.5. CO 

We observed two settings in the 5 /im v — 1-0 band of 
CO, which included a number of P-branch lines of sev- 
eral isotopologues. These settings were centered approx- 
imately at 2055 cm~^ and 2085 cm~^. Strongest among 
the detected features were the ^^CO lines with J = 14, 
15, 21 and 22. Each line exhibited a strong P Cygni 
profile with emission and absorption components sepa- 
rated by ^20 km s~^. Each component was broader than 
those of any line of any other molecular species we ob- 
served, with line widths of 10-20 km s~^. In each case the 
strength of the emission was approximately equal to the 
strength of the absorption. The features are consistent 
with an expansion or outfl ow and have been previously 
interpreted as such by Mitchell & Hasegawa (1991). Ob- 
servations of IRS 1 were also made at these settings; we 
find a different profile altogether with a double-bottomed 
absorption feature in ^^CO, weaker absorption on the 
blue wing and a slight emission bump on the red wing 
of each feature. Spectra of both objects are shown in 
Figure [7[ 

Four^^CO P-branch lines with J = 3, 4, 11, and 
12 were also detected in the spectrum of IRS 9. Of 
these, P(ll) was a clear detection near the middle of 
a spectTral order and away from contamination by tel- 
luric atmospheric lines, P(12) was noisy and found very 
near the edge of a spectral order, and P(8) and P(9) are 
clearly detected, but blended slightly with two nearby 
C^^O features. Additional "CO lines, P(ll) and P(12), 
are located in the spectral setting at 2055 cm~^ but 
occur in regions of strong telluric atmospheric absorp- 
tion and thus were not detected. We also searched for 
lines attributable to C^^O but did not confidently detect 
any; this is not surprising giv en the low expected abun- 



Smith et al. 



(2009) quote 



dance of C^^O. For reference, 

jY(12q16q-)/^(12(2;17o) of 2800±300 towardlheT Tauri 
star VV CrA, obtained at a spectral resolution similar 
to our TEXES observations. They compare this to the 
Local ISM value of 2005±155 ( |Wilson|p;999| ) 

The CO lines seen toward both IRH 9 and IRS 1 have 
different line shapes in ^^CO and ^-^CO. In the case of 
IRS 1, the ^^CO features show only single components 
in the deepest absorption, while toward IRS 9 the ^^^CO 
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Fig. 6.— TEXES spectra toward NGC 7538 IRS 9 at the position 
of the ^^NHs U2 antisymmetric Q branch (top panel) and symmet- 
ric P{7,K) branch (bottom panel). The upper curve in each plot 
shows the data, while the lower curve in each shows the relative 
atmospheric transmission on the same scale. Asterisks mark posi- 
tions of the NH3 lines. Both spectra have been corrected for the 
Earth's motion relative to the LSR, and the position of the NH3 
line markers shifted to the systemic LSR velocity of approximately 
-60 km s~^. 

P(ll) line and ^^CO lines show the pronounced P Cygni 
profiles. The ^'^CO lines in the 2085 cm^^ setting clearly 
show multiple components suggestive of P Cygni pro- 
files in which emission is weak relative to absorption, 
and there are at least two distinct absorption units. In 
turn, the strongest absorption components show a weak, 
double-bottomed structure like the ^•^CO lines in IRS I. 
The multiple absorption units may appear in the profiles 
of the ^^CO features in the same spectral setting; weak 
"shoulders" on the blue wing of the absorption compo- 
nent are consistent with additional absorption suffering 
contamination by the primary components. Such shoul- 
ders are also noticeable on the red wings of the emission 
features of ^^CO . These details can be seen in Figure [sj 

The detection of ^-^CO toward IRS 9 with line depths 
comparable to those of ^^CO strongly indicates satura- 
tion of both isotopologues. Useful columns may how- 
ever be obtained from rarer variants such as C^^O, but 
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Fig. 7. — The upper panel shows spectra of two P-branch lines of 
^^CO and one P-branch line of ^^CO toward IRS 9 obtained with 
TEXES. The lower panel shows the same spectral settings observed 
toward IRS 1. In both cases the motion of the Earth with respect 
to the LSR has been removed from the velocity scales. 

these lines are too weak to be measured accurately in 
the present spectra. It would also be of interest measure 
the ^^C/^^C ratio toward IRS 9 for c omparison w ith the 
value of 45 toward IRS 1 reported by |Knez et al.| p009). 
As we did not observe the ^^C isotopologue of C2II2 in 
our spectra of IRS 9, the CO observations would offer a 
way to measure the ratio if the lines were not saturated. 

3.6. OCS 

The ^^CO setting at 2055 cm~^ we observed also con- 
tains the positions of a number of P-branch lines of the 
1^2 stretching mode of OCS (12 < J < 22), but we can- 
not confidently claim detections of any of them toward 
IRS 9. Periodic structure in the continuum of our 5 /im 
spectra may be broad OCS features, although it is more 
likely the structure results from flat- fielding residuals. 
For reference, we have overplotted the locations of the 
rest-frame OCS line centers over the data in Figure [9 
It is unclear whether we should have observed lines of 
OCS in the gas phase, given the strong detections of 
CO nearby. Low temperatures may inhibit desorption 
of OCS from grains, but the theoretical gro unds for pre- 
dicting OCS abundances are questionable (Doty et al 
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Fig. 8. — TEXES spectrum toward NGC 7538 IRS 9 centered near the 5 fira rotational lines of several isotopologues of CO. The data are 
shown in the upper trace (solid line) and the relative atmospheric transmission in the lower trace (dot-dashed line) multiplied by a factor 
of 10. Some line identifications are shown. 

was devised in which rough estimates of the parameters 
were made from "by eye" fits of synthetic spectra to our 
data; these values were then used as the initial guesses 
in our fitting program which iterates values of hundreds 
parameters to achieve the best fit to the data by minimiz- 
ing x^- While our code is robust in finding solutions 
even when starting far from the global minimum, con- 
vergence is achieved faster by using good initial guesses. 
This approach was best suited to our method of spectral 
analysis. 

4.1. Simple Synthetic Spectrum Fitting 

In order to get approximate parameters for the gas- 
phase lines we observed, we first fit synthetic spectra to 
the data using a fitting method commonly referred to 
as "chi by eye." This approach allowed us to get quick 
initial estimates of parameters such as excitation temper- 
ature and column density for the various molecules that 
were used as inputs to a more sophisticated non-linear 
least squares code to obtain optimal values of the param- 
eters. The simple method also allowed us to determine 
the number, width and velocities of line components. 

The equivalent width of an optically-thin line is the 
product of the line strength, a, and the column density, 
N. We obtained values of a for each molecular species 

of laboratory molec- 
2005;). The GEISA 



20041. The signature of an outflow in the spectra of CO 
hints at possible excitation through shocks. We might 
expect that OCS would be similarly excited but there is 
no support f or this in our data . This is contrasted with 
the results of Evans et al. ( 1991), who detected gas-phase 
OCS toward Orion iKc2; for reasons given in Section 5, it 
is impossible to know whether this is due to temperature 
differences or real abundance variations between the IRS 
9 and IRS 1. 

3.7. HNCO 

We looked for 1/4 bending mode HNCO IL-0 P-branch 
lines toward IRS 9 in the TEXES spectral setting that 
included s(7,A') lines in the P branch of ammonia. None 
of the three lines included in this setting (2 < J < 4) 
was detected, although there may be HNCO Q-branch 
absorption at 827.7 cm~^, obscured by telluric interfer- 
ence. Given the strong observed ice ab sorption at 4.62 
/xm toward IRS 9 ( Pendleton et al.|1999 ) , now commonly 
attributed to OCIN", and a lack or gas-jDhase HNCO ab- 
sorption, it seems most likely that temperatures in the 
outer envelope of IRS 9 are sufficiently low to prevent 
OCN~ desorbing from grains in detectable quantities. 
Alternately, gas-phase OCN~ might exist at H (or H2) 
abundances insufficient to result in considerable HNCO 
production. 

4. ANALYSIS AND MODELING 

We fit synthetic spectra to the TEXES data to de- 
rive excitation parameters and column densities for the 
various molecular species observed. A two-step process 



observed fr om the GEISA databas^ 



ular data ( |Jacquinet-Husson et al 



coefficients are specified tor a fiduciaTtemperature — 
296 K; we corrected them to an assumed temperature T 

^ http: / /ether. ipsl.jussiou.fr/ethcr/pubipsl/geisa_iasi_raie_frame_2003_uk.jsp 
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Fig. 9. — The TEXES spectra of NGC 7538 IRS 1 (top trace in each panel) and IRS 9 (bottom trace) near features in the P-branches 
of i^CO and ^^CO; heavy solid lines are the object spectra and light dot-dash lines are the sky spectra offset for clarity. The wavenumber 
scale has been corrected to a velocity reference at rest in each object's frame, and the line centers of OCS features, marked with asterisks, 
have been velocity shifted appropriately. 
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via the formula 



a{T) — a{To) exp 



(2) 



where Ei is the energy of the lower state of a transi- 
tion, ks is Boltzmann's constant, and /3 = 1 for linear 
molecules and | for non-linear molecules. Here the tem- 
perature dependence of the vibrational partition function 
and the stimulated emission correction are neglected as 
both are of order unity. For Gaussian lines, the line- 
center optical depth is calculated according to 



n 



Na 



(3) 



where b is the Doppler b parameter in wavenumber units. 
The run of optical depth with frequency over the line is 
computed with 



t{w) = Ti exp ■ 



-(w - Wq)^ 



(4) 



in which w is the wavenumber and Wo is the line-center 
wavenumber. The optical depths are converted to in- 
tensities by taking the exponential of t{w); the contin- 
uua are normalized to unity. We assume that the ro- 
tational states are in local thermodynamic equilibrium 
(LTE). Our radiative transfer model also assumes pure 
absorption for fitting the observed absorption lines. Re- 
emission by the absorbing molecules is neglected. Emis- 
sion features are treated as separate components, which 
are modeled as negative absorption. Populations of ex- 
cited vibrational states are not calculated. Temperatures 
are arbitrarily held fixed at 200 K; this value is consistent 
with temperatures around IRS 1 derived from TEXES 
data. The exception was the Q-branch of C2H2 where a 
number of observed J states allowed the temperature to 
be constrained to 100 K. Column densities and velocities 
derived using this method are shown in Table [2j Fits to 
the individual lines are shown in Figures [T0| through [12) 



4.2. Minimization Spectral Fitting 

Our spectral fitting code, FITSPEC, uses the Levenberg- 
Marquardt method of inv erting the second derivative ma- 



trix to minimize a-s in Bevington & Robinson ( 2003 1 
for up to four unique absorbing or emitting units, or 
"clouds," along the line of sight. The 1-D radiative trans- 
fer model assumes rotational LTE of all molecules with a 
single excitation temperature. As in our "by eye" fitting, 
re-emission by absorbing molecules is neglected and emis- 
sion features are modeled as negative absorption. The 
total number of free fitting parameters depends on the 
number of molecules in the fit and the number of ech- 
elon orders in a particular spectral setting. Although 
there is a large number of free parameters in the fit, most 
of them are "uninteresting," such as baseline coefficients. 
The "interesting" parameters in the fit are Doppler shift, 
Doppler b parameter, column density, temperature, and 
covering factor for each cloud and for each molecu l e. Th e 
details of the program are discussed in Knez et al. ( 2009 ) . 
For our observed molecules, we chose a number of absorb- 
ing clouds consistent with the number of velocity compo- 
nents seen in the TEXES spectra, allowing FITSPEC to 



independently determine their parameters. When avail- 
able, we included data indicating the non-detection of 
certain lines, such as high-J C2H2 near 13 fim, to better 
constrain the fit parameters. 

The success of FITSPEC requires accurate molecular 
transition data to compute synthetic spectra against 
which to fit the d ata. We relied on the calculations in 
Evans et al.] ( |1991[ ), giving information including the line 
center, the energy of the transition either in units of J( J 
-1-1) (which is then multiplied by the molecule's rota- 
tional constant B) or E in units of cm~^, and the ro- 
tational line strengths (Honl-London factors) in units of 
cm~^ (10^^ cm"^)"-'^. The program uses guesses of the 
starting parameter values to compute the initial model 
spectrum, and then iterates the fit by varying those val- 
ues until a convergence criterion is satisfied. The require- 
ment is 



< 0.01 



(5) 



where Api is the change in the value of parameter i and 
(Tp. is the uncertainty in the parameter i, allowing all 
other parameters to vary. 

The output of the program consists of the best-fit pa- 
rameter values, uncertainty estimates and a synthetic 
spectrum having the properties of the best fit. The in- 
tensity of the synthetic spectrum is calculated by taking 
I = Ioe~'^ of each absorbing unit along the line of sight 
and adding any optically thin emission as 

N 

lobs = /O n [1 - ^"(1 - e"""^"^)] + SemCemTemH (6) 
n=l 

where N is the number of absorbing clouds in the fit, 
Iq is the intensity of radiation from the background con- 
tinuum source (set to unity since the spectra are nor- 
malized), Cn = and r„ are the covering factors 
and optical depths of the absorbing clouds, Sem and 



Ce 



are the source function and covering fac- 



tor of the emitting material, and Hq is the solid angle of 
the continuum source. The right-hand term represents 
an optically-thin emitting cloud; since it only depends 
on the product of Sem and Cem, we scale Cem by Sem/Io 
and drop Sem from the formula. The product term in 
Equation [6] assumes that the absorbing clouds overlap 
randomly. That is, if two clouds each have covering fac- 
tors of 0.1, 1/10 of each lies along the line of sight to the 
other. 



Knez et al. (20091 used the same approach and note 
this results in greater absorption in saturated lines than 
would result if the optical depths were added first and 
then the transmission spectrum were calculated from the 
optical depth spectrum. The resulting fits to the data 
are better than those obtained by summing the optical 
depths first and then taking the exponential. This fitting 
approach, however, is not without drawbacks. Emission 
lines may not be thermally excited, but rather could be 
due to resonantly scattered light; we have no way of dis- 
criminating between these processes. We assume that 
upward transitions in the molecules we observe are not 
followed by re-emission of that radiation; in reality, most 
absorptions should be followed by re-emission, although 
in a direction uncorrelated with the incoming radiation. 
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TABLE 2 

Molecular absorption parameters derived from simple synthetic spectrum fitting in LTE 



JVIoIscuIg/ rjcLriQ 




a 

n 


T^ 


f 1 

(km s^^) 


f 2 


^3 


(km s~^) 




"3 


iVi 

(cm^2^ 


iV2 


iV3 


C2H2 1^5 


R{1) 


1 


100 


-62.0 






3.3 






1.0(15) 








R{5) 


1 


100 


-62.0 






2.4 






4.5(14) 








R(6) 


1 


100 


-65.5 






2.4 






4.5(14) 








Q branch 


1 


100 








1.0 






2.0(15) 






C2H2 U4, + 1/5 


P(2) 


1 


100 


-61.0 






6.8 






2.5(16) 








P(3) 


1 


100 


-62.0 






5.7 






1.0(16) 






HCN 1/2 


R(16) 


1 


200 


-64.5 


_ 




7.9 


_ 




1.2(16) 


_ 




12CH4 1/4, 


R(O)'^ 


4 


200 


-65.0 


-54.0 


-61.0 


2.3 


2.7 


4.6 


7.5(16) 


-4.0(15) 


6.0(16) 




R(2) 


2 


200 


-65.5 


-61.5 




2.3 


2.5 




6.0(16) 


3.3(16) 






R{2) 












2.5 


2.3 




6.3(16) 


1.5(16) 




"NH3 U2 


aP(4,0) 


2 


200 


-64.0 


-70.0 




4.6 


3.9 




-2.5(15) 


5.0(14) 






aP(4,l) 


2 


200 


-62.0 


-72.0 




5.3 


3.9 




-3.8(15) 


2.5(14) 






aP(4.3) 


2 


200 


-62.0 


-70.0 




2.8 


3.9 




-1.0(15) 


1.0(15) 






sP(7,0) 


2 


200 


-60.0 


-69.0 


_ 


2.9 


2.9 




-1.5(15) 


1.2(15) 






sP(7,l) 


2 


200 


-62.0 


-70.0 


_ 


1.8 


4.0 




-2.0(15) 


7.5(14) 






sP(7,2) 


2 


200 


-61.0 


-62.5 


_ 


3.3 


2.9 




-6.0(15) 


4.5(15) 






sP(7,3) 


2 


200 


-58.5 


-64.0 


_ 


1.8 


2.9 




-1.2(15) 


1.5(15) 






sP(7,4) 


2 


200 


-60.5 


-63.0 


_ 


1.8 


1.8 




-2.5(15) 


2.0(15) 






sP(7,5) 


1 


200 


-56.0 


- 




4.4 


- 


- 


-7.0(15) 


- 


- 




sP(7,6) 


1 


200 


-57.0 






4.4 






-4.0(15) 






12C160 v=l-0 


P(14) 


3 


200 


-55.0 


-64.5 


-88.5 


2.4 


6.5 


7.2 


-1.0(16) 


-1.1(17) 


3.0(17) 




P(15) 


3 


200 


-51.5 


-63.5 


-88.5 


3.6 


9.4 


7.2 


-2.5(16) 


-2.2(17) 


3.0(17) 




P(21) 


2 


200 


-54.5 


-77.5 




10.9 


10.9 




-8.0(17) 


5.5(17) 






P(22) 


2 


200 


-55.0 


-77.0 




8.0 


10.2 




-6.0(17) 


6.0(17) 




iSCif^O v=l-0 


P(3) 


3 


200 


-60.5 


-70.5 


-89.0 


4.3 


4.3 


1.4 


-9.0(17) 


8.0(18) 


7.0(17) 




P(4) 


3 


200 


-61.5 


-71.5 


-86.0 


3.6 


5.0 


6.5 


-7.0(17) 


6.8(18) 


3.4(18) 




P(ll) 


2 


200 


-54.0 


-75.0 




5.8 


10.2 




-8.0(17) 


2.7(18) 






P(12) 


2 


200 


-54.0 


-73.0 




5.9 


10.2 




-3.5(17) 


6.5(18) 




12C180 v=l-0 


P(2) 


2 


200 


-68.0 


-72.0 




2.9 


3.6 




7.5(18) 


7.5(18) 






P(3) 


2 


200 


-66.0 


-71.0 




4.3 


4.3 




3.0(18) 


3.0(18) 





Number of velocity components fit. 
^ All velocities are referred to the LSR. 

Column densities are given in the form A{B) = Ax 10^. Negative values of A indicate components in emission. "Column density" 
includes all rotational states, assuming T = 200 K. 

The best fit to this line required four components. The remaining component not in the table is V4 
s-i, = 5.0(15)cm-2 



-71.5 km s-i, bi = 2.3 km 



SO the extent to which re-emission fihs in the absorption 
depends on the distribution of the re-emitting gas around 
the continuum source. 

The results of our FITSPEC runs on the IRS 9 data are 
summarized in Table |3] where they are compared with 
results from fitting TEX ES data of IRS 1 with the same 



where x = v — v' . The upper limit of the equivalent 
width, Wmax, of a line not detected is proportional to 
the peak value of the convolved spectrum according to 



W„ 



(9) 



procedure in Knez et al. (2009). 



4.3. Column Density Upper Limits For N on- Detections 

We calculate column density upper limits from the 
lines that were not detected in the spectrum of IRS 9 by 
convolving the data with a Gaussian line shape function 
whose properties are consistent with those of detected 
lines. We begin by assuming a normalized spectrum con- 
taining a single absorption line has the functional form 



(7) 



where a > is the depth of the line and h is the Doppler 

h parameter. The observed spectrum is convolved with a 

2 /1 2 

Gaussian of the form ne~^ ' where n is a normalization 
factor here chosen to be 1/ ^Jlxh such that the integral of 
the Gaussian is unity. The convolution is given by 



> ne 



^y"'- = / F^ne-^'/^'d., 



(8) 



in which P is the peak. We used 6 = 3.5 km s^^ in making 
the calculations, a value consistent with those derived 
by spectral fitting methods described previously. The 
column density limit for species X was then computed 
from 

N(X) = (10) 
a 

where a is the line strength. We used the GEISA line 
strengths where available and corrected them to a uni- 
form temperature of 200 K as described previously. Re- 
sults for our non-detections are given in Table |4j 

We derived upper limits for all lines of a particular 
molecular species available in a given spectral setting, as 
indicated in the table by the number of lines parameter 
n. In cases where n > 1, the equivalent width and column 
density upper limits given are simple averages. All limits 
are given with respect to the total number of molecules 
along the column, irrespective of isotopic composition 
except for ^'^CH4; here we have divided the derived col- 
umn density by the terrestrial ^"^C/^^C = 89 ratio built 
into GEISA to recover the column density specifically for 
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Fig. 10. — T EXES spectra of NGC 7538 IRS 9 (thin lines) and synthetic spectra computed using the simple LTE fitting method described 
in Section |4.1| (thick lines). The continuum in each spectrum has been normalized to unity and the spectra are plotted on a velocity scale 
referred toTTie LSR. Features shown are C2H2 U5 R{5), R{6) and iJ(13); C2H2 U4 + 1/5 P(2) and P(3); HCN U2 i?(16); and CH4 U4, R{0) 
and the two components of R(2). 
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4.4. An Upper Limit to the Ionizing Flux in IRS 9 

The ^P3/2 fine structure line of [Ne II] at 

780.42 cm^^ is often used as a diagnostic of ionized gas 
at our wavelengths. The non-detection of this line in 
the spectrum of IRS 9 places a constraint on the amount 
of ionizing flux being emitted by the embedded central 
source. 

The Lyman continuum luminosity (N^yc) can be cal- 
culated from the [Ne II] luminosity since both are pro- 
portional to the square of the density integrated over 
volume. The [Nc II] volume emissivity is given by 



Tf, = 10"' K) is the critical densi ty in the two-level ap- 
proxim ation, using the values in |Osterbrock fc Ferland| 
( 2006 1 . The hydrogen recombination rate per unit vol- 

where as (= 



ume IS given by n. 



xlO" 



-13 



recomb — O^QTl^Tlp 



2.6 



cm'^ s~^) is the recombination rate for "case B" 
recombination, and rip is the proton density. The ratio 
of N^yc to the [Ne II] luminosity is then given by 



Lyc 



^recomb 



nH+ OLB 



L[NeII] j[NeII] «Afe+ qiuhl' 



(12) 



j[Neii] = n^nN^+qiuhv- 



(11) 



where rig is the electron density, n^ve-*- is the density of 
ionized neon, qi^ is the rate coefficient for coUisional ex- 
citation {qiuhv — 8.4 X 10~^^ erg s~^ cm^), and = 
Aui/qui (= 4.9 X 10^ cm~'^ at an electron temperature 



assuming ^ ric- The Ne"'"/H+ ratio is given by the 
product of the Ne abundance and the singly-ionized frac- 
tion for Ne, ,f+. We ta k e the value of /+ — 0.3 from 
Figure 3 in Ho & Keto (20071 using for the radiation 
temperature that of a /AlvlH star with IRS 9's luminos- 
ity {Tr — 26200 K) and a ratio of the dilution factor for 
blackbody radiation, F, to the electron density of F/n^ 
= 10^^^ cm'^. From this we adopt a Ne+ / 11+ ratio to 
be 3 X 10"^ 

Without knowing the width of the [Ne II] line, it is 
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Fig. 11.— TEXES spectra of NGC 7538 IRS 9 (continued). Features shown are NH3 1^2 aP(4,0), aP(4,l), aP(4,3), sP{7,0), sP{7,l), 
sP(7,2), sP(7,3), sP(7,4)i sP(7,5) and sP(7,6). The sP(7,0) and sP(7,l) lines are shown together in a single panel using the velocity scale 
for sP(7,0), and are offset slightly for clarity. 



difficult to place a strong upper limit on its flux, but 
an emission equivalent width of Wi, = 0.01 cm^^ 
be ruled out unless the line is wider than 0.2 cm 
or 77 km s~^. This corresponds to a line flux limit of 



can 

1 



Nell] 
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^ and a luminosity limit 
IQ'i^ era; s We then 



< 2.0 X 10 erg s " cm 
of i[Are//] = ATTd'^F[Neii] < 1-7 X lU"' erg s 
calculate NLyc < [3.3 x 10"* x 2.6 x 10-1^8-4 x 10" 
1.7 X 10^^ = 1.7 X lO'*^ photons s^^. For reference. 
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lens (2005) gives the Lyman continuum luminosity of a 
main sequence BO star as N^yc = 1.4 x 10"*^ photons s~^. 
Apparently IRS 9 emits a small fraction of the Lyman 
continuum luminosity of a main-sequence star of its bolo- 
metric luminosity. The radiation environment of IRS 1 is 
considerably more intense, by comparison. From infrared 
and millimeter data, the exciting source is inferred to be 
of spectral type earlier then 07.5, total luminosity L > 
8 X 10^ Lq, and Lyman continuum luminosity iV/,,,^ > 



lO''^ s-1 (Hackwellet al.|l982 


Akabane et al.|2001 


Lugo 


et al.|2004p. I'he interred radiation held ot IRS 1 is more 



of conditions under which chemistry proceeds. 



5. INTERPRETATION OF THE SPECTROSCOPIC RESULTS 

5.1. Assumptions and Caveats 

Interpretation of our results requires caution, given the 
underlying assumptions of the radiative transfer model 
in FITSPEC. The program assumes a single continuum 
source behind all absorbing material and a constant tem- 
perature characterizing the material in each of up to four 
absorbing "clouds" along the line of sight. We also add 
an optically-thin emitting cloud, but neglect re-emission 
by the absorbing gas and absorption by the emitting gas. 
The column densities we report do not represent the en- 
tire column down to the embedded luminosity source, 
but rather the column only to an effective "dust photo- 
sphere" interior to which the continuum optical depth is 
large at a given wavelength. Consequently, our observa- 
tions at 5 fim and 8 fim likely probe to greater depth 
than those at 13 /im. The optical depth may be higher 
at longer wavelengths due to silicate opacity, and the 
colder, outer layers of the dusty region are likely too 
cool to emit significantly at 5 /im. Comparison of results 
at different wavelengths is therefore difficult. We also 
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Fig. 12.— TEXES spectra of NGC 7538 IRS 9 (continued). Features shown are ^^C^^O v=l-0 P(14), P(15), P(21) and P(22); i^ciBo 
^=1-0 P(3) / i^cISq v=1-0 P(2); ^^C^^O d=1-0 P(4) / ^^C^^O v=l-0 P(3); and "cISq i'=1-0 P(11) and P(12). 



did not take into account emission by molecules follow- 
ing absorption of photons, except as approximated by 
our "negative absorption" method of treating emission 
lines. Furthermore, we ignored emission following radia- 
tive excitation. We assumed either pure absorption or 
optically-thin emission. However, the neglect of emission 
from a radiatively excited molecule should not be valid; 
at any plausible gas density the coUisional de-excitation 
rate should be far below the radiative rate. Critical den- 
sities for vibrational transitions are typically greater than 
10^^ cm~'^, and almost all absorption events should be 
followed by emission, not collisional de-excitation. For 
two of our molecules, NH3 and C2H2, we were able to 
observe lines of more than one vibrational band or rota- 
tional branch, and the comparison of the observed line 
depths provides evidence that re-emission following ab- 
sorption is an important effect. 

For NH3 , we observed weak emission in P-branch lines 
at wavelengths near 12 fim and possibly weak absorp- 
tion in Q- branch lines near 10.5 /im. i?-branch lines 
were not observed due to the silicate and telluric O3 ab- 
sorption at their wavelengths. We explain the differences 



between the P- and Q-branch lines based on the differing 
continuum radiation fields due to their differing proxim- 
ity to the silicate dust feature. The NH3 R branch is 
centered on the 9.7 /im silicate feature, the Q branch 
falls on its long wavelength shoulder, and the P branch 
is mostly between the 9.7 and 18 /im features. If the 
radiation field seen by a molecule is due to silicate emis- 
sion cither because the molecule is exposed to radiation 
from an optically thin external dust cloud or because it 
is immersed in a cloud of warm silicate dust, there will 
be more upward radiative transitions (i.e. more absorp- 
tion) in the Q and R branches than in the P branch. 
The branching ratio between P-, Q-, and i?-branch tran- 
sitions from a given vibrationally-excited state depends 
on the K rotational quantum number, but is similar in 
the three bands. Consequently, if all absorption events 
are followed by an emission, absorptions will outnumber 
emissions in the R branch and emissions will outnumbers 
absorptions in the P branch, resulting in net emi ssion in 
the P-branch lines, as we observe. In Section |5.3| we 
describe a simple model involving purely radiative tran- 
sitions in NH3 incorporating upward and downward line 
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TABLE 3 

Molecular absorption parameters determined from 
minimization spectral fitting 



Species vj^sR "^ex 
(km s-i) (K) 



b N(xy 
(km s-i) (10^® cm 



C 

% 



IRS 9 (this work) 



C2H2 


-58.9(1.4) 


195(57) 


4.0(1.3) 


-0.2(0.1) 


52(189)*^ 




-63.9(1.0) 


205(49) 


4.0(1.1) 


2.3(0.7) 


16(9) 


HCN 


-58.7(1.7) 


200(68) 


4.0(1.3) 


0.5(0.7) 


52(193)'' 




-63.8(2.9) 


204(76) 


4.0(1.4) 


3.9(4.3) 


11(6) 




-57.3(0.5) 


365(47) 


5.1(0.7) 


-1.1(0.1) 


52(195)'' 




-61.6(0.6) 


201(43) 


5.9(1.1) 


4.8(0.8) 


8(8) 


HNCC^ 


-58.9(0.0) 


195(0) 


4.0(0.1) 


0.1(-0.5) 


52(0) 




-63.9(10.0) 


205(50) 


4.0(1.0) 


0.9(-1.6) 


16(10) 



IRS 1 (Kncz et 



al.|2009[ l 



C2H2 


-55.7(0.3) 


225(20) 


0.6(0.1) 


3.0(0.6) 


6(0.6) 




-59.4(0.3) 


191(10) 


0.6(0.1) 


2.8(0.2) 


24(1) 


HCN 


-56.3(0.3) 


256(30) 


0.8(0.1) 


5.6(2.8) 


8(1) 




-60.0(0.3) 


456(127) 


1.0(0.4) 


1.3(1.2) 


68(50) 


NH3 


-57.3(0.3) 


278(31) 


0.9(1.0) 


5.2(1.5) 


9(1) 




-60.1(0.3) 


248(32) 


0.4(0.1) 


2.8(0.8) 


16(4) 


HNCO 


-57.2(0.3) 


319(27) 


1.8(0.2) 


0.4(0.3) 


4(0.3) 




-60.2(0.2) 


171(79) 


0.7(0.4) 


0.1(0.3) 


97(11) 



Note. — Quantities in parentheses are the la uncertain- 
ties on the adjacent figures in the same units unless otherwise 
indicated. 

^ Negative column densities indicate components in emission. 
Column densities for the emission components are Cem''"em('^) 
converted to and do not represent a real column densities. 
^ Unable to be constrained in the fit, assuming lines are opti- 
cally thin. 

All HNCO parameters other than the column density were 
fixed equal to those of C2H2 because an insufficient number of 
HNCO lines were observed to constrain them. 



TABLE 4 

3a UPPER LIMITS TO COLUMN DENSITY OF MOLECULES IN 
THE UPPER STATE FOR SELECTED NON-DETECTIONS 



Species/Mode 


Line 




NumitiX) (cm-2) 


"CH4 U4^ 


i?(l) 


1 


4.1x101^^ 




R{4) 


4 


3.1xl0i« 


OCS U2 


P branch 


34 


2.9x1015 


HCN U2 


R{15) 


1 


1.8x10^'' 




R(22) 


1 


8.5x10^'' 




R(23) 


1 


1.1x10^^ 


C2H2 1^5 


R{12) 


1 


8.5x10^5 




R(21) 


1 


9.4x1015 




R(22) 


1 


3.5x10^^ 


C2H2 21/5 - 1/5 


Q branch 


15 


2.1x10^^ 



Note. — A uniform temperature of 200 K was assumed 
for all species. 

^ Number of lines/components in the attempted fit. 

The given values of A(i^CH4) are for the isotopologue 
only rather than the total CH4 column density. 

strengths and assuming that the incident radiation field 
is proportional to the strength of the silicate feature at 
the wavelengths of the various lines. The model is time- 
independent and assumes the NH3 gas is optically thin. 
More detailed and quantitative models will follow in a 
future paper. 

5.2. C2H2, HCN, HNCO, OCS, and CO 



We find that excitation temperatures for various 
molecules are basically consistent toward IRS 9 and IRS 
1 with the exception of the hotter of the two NH3 compo- 
nents toward IRS 1. That value of T^x has a rather large 
uncertainty so it may in fact lie closer to ^250 K. The col- 
umn densities are also similar in the two sources, though 
again with p articular exception in the case of HCN. [Knez| 



et al. (2009) had the benefit of more observed HCN lines 
but we find a lower column toward IRS 9. Most striking 
is the disagreement between the two sources in terms of 
the intrinsic line widths; although the thermal broaden- 
ing should be similar, turbulence or systematic velocity 
variations along the line of sight to IRS 9 are greater 
than toward IRS 1. At the least, we can say that the 
spectral lines observed in the direction of IRS 9 are gen- 
erally less saturated than those seen toward IRS 1, and 
this has important consequences for the accuracy of the 
IRS 9 column densities we obtain. 

An example in which saturat ion of lines becomes im- 
portant is CO. Knez et al. ( 2009 ) derived fractional abun- 
dances relative to CO for a variety of molecules seen 
in absorption toward IRS 1, attempting to compensate 
for badly s aturated ^^CO lines by using the value of 
N(}^CO) in |Mitchell eral] (|l990') and adopting ^"^C/^^C 
= 45 on the basis of their observations of C2H2. The 
value of N{^^CO) for IRS 9 reported by MitcheU et al. 
was based on their measurement of ^^C O lines and as- 
sumed the terrestrial ratio -'^C/"''^C = 89. Boogert et al. 
(20041 confirm the saturation of the ^^CO lines observed 



toward IRS 9. From observations of lines of ^'^CO and 
12q18q upper J levels of 7-15 they derive an opti- 

12C/13C = 80 from 



cal depth T^2 = 13 ib 4 assuming 



Boogert et al. ( |2002 ). However, they argue that reliable 
valuesforlVp^CO) can be obtained, since the ^'^CO lines 
remain optically thin; they quote a value A^("'^^CO) = 
(3.2±1.0)x lO^'* cm^2. While we may argue with their 
choice of ^^C/^^C and with their conclusion that the 
^■^CO lines are optically thin, their value of the CO col- 
umn density is probably more realistic for IRS 9 than 
that obtained from t he lower-resolution observations of 
Mitchell et al. (1990). Abundance ratios for some of the 
molecules we observed with TEXES, relative to this col- 
umn of CO, are shown in Table [5] 

Similarly, an estimate of the abundance of various 
molecules with respect to H2 can be made given a mea- 
surement of the quantity of dust along the line of sight 
and a value for the gas-to-dust ratio in a typical proto- 
stellar envelope. Knez et al. adopted an H2 column den- 
sity toward IRS 1 of 7.5x10 ^^ cm"^ based on the 9.7 ^lui 
optical depth measured by Willner et al. (1982) to cal- 
culate fractional abundances of other molecules. Using 
the same method, we derive 7V(H2) = 5.2 x 10^^ cm"^ 
toward IRS 9 and present the resulting abundances in Ta- 
ble [5j The pattern of relative abundances differs between 
the objects according to both the molecular species and 
whether the abundance of a given species is referred to 
CO or to H2. The latter observation indicates that either 
the abundance ratio of CO to H2 or the gas-to-dust ratio 
varies between objects. 

The assumptions made in calculating the values in 
Table |5] - that the CO column density is comparable 
along both lines of sight and that the dust distribu- 
tion and gas-to-dust ratio are similar in both cases - 
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TABLE 5 

Abundances in IRS 9 and IRS 1 with respect to CO and H2 



Molecule 


IRS 9 (this work) 


IRS 1 (|Knez et al.|2009|l 


N{X)/N{CO) 


iV(X)/Ar(H2) 


Ar(X)/JV(CO) 


iV(X)/iV(H2) 


C2H2 


7.7(-3) 


3.3(-7) 


5.7(-3) 


7.7(-7) 


HCN 


1.4(-2) 


5.8(-7) 


6.9(-3) 


9.2(-7) 


NH3 


1.9(-2) 


7.9(-7) 


3.6(-2) 


4.8(-6) 


CH4 


9.9(-2) 


4.2(-6) 


3.6(-2) 


4.8(-6) 


HNCO 


<3.0(-3) 


<1.3(-7) 


5.0(-4) 


6.0(-6) 



Note. — All values are given in the notation A{B) = ^ X 10^. The c olumn 
densities to ward IRS 1 adopted for CO and H2 were Ar(i^CO) = IxlQis cm-2 | |Knez| 
|et al.|2009| and N(}i2) = 7.5x10^2 ^.^-2 jwillner et a l.|l982t. For IRS 9 t he CO 
column density used was N{^^CO) = 3.2x10^° c m~^ (Boogcrt et al.|2002| and the 
H2 column density was the same as that used by |Knez et al.| ( |2UU9( I . 



may not be unreasonable given that both objects are 
along nearby sightlines toward the same molecular cloud. 
Given this caveat, it is worth noting that the fractional 
abundances of all molecules relative to both CO and H2 
are broadly consistent with each other, except for HNCO 
with the values given for IRS 9 properly unconstrained by 
FITSPEC because no HNCO hues were detected. C2H2, 
HCN, and NH3 are generally less abundant in IRS 9 than 
IRS 1, although this is only a suggestion in light of the 
assumptions. 

Some kinematic information about IRS 9 is available in 
the line shapes. In combination with the varying effective 
depths to which sightlines probe at different wavelengths, 
line shapes give a sense of the distribution of various 
molecules along the line of sight. The P Cygni line pro- 
files of our CO observations, for example, clearly indicate 
an outflow known previously. CH4 absorption toward 
IRS 9 was known previously as well, but the TEXES 
data reveal the structure of the lines in more detail than 



the Keck+NIRSPEC spectr a of|Boogert et ar](|2004| ) and 
the IRTF+IrsheU spectra of |La!cy et al.| ( |199 1|). 



In contrast to the circumstellar disk model of IRS 1 



presented in Knez et al. (20091 based on TEXES obser- 
vations, we cannot conciusiveiy establish the presence of 
a disk in IRS 9. However, our fitting program preferred 
relatively small covering factors, similar to those in the 
best-fit synthetic spectra for the lower velocity compo- 
nent of IRS 1. Knez et al. interpret this as the signature 
of a near edge-on disk, acting as a source of continuum 
radiation that is only partially covered by the absorbing 
material. Given the high- velocity out flow in IRS 9 re- 
ported by Mitchell & Hasegawa (19911), we might rather 



be observing an accreting system nearly face-on in which 
the outflow both sweeps out a region near the contin- 
uum source and entrains material into the flow that only 
partially covers the source. Thi s picture would accou nt 



for the methane observations of Boogert et al. ( 2004 ) in 
which cold, solid-phase methane is seen in the outer en- 
velope, becomes depleted in the outflow, and is seen in 
the gas phase in the inner envelope, all along a s ingl e 
sightline. We develop this model further in Section [5.6[ 
Non-detections of molecules also furnish some inter- 
pretable information. Carbonyl sulfide (OCS) is a 
molecule that commonly occurs in ices along sightlines 
through molecular clouds and indicates the presence of 
sulfur-bearing precursor species such as H2S and SO2 



( [Ferrante et"aI]|2008L |Gibb et all (l2004[ ) detected sohd- 
phase OCS toward IRS 9 with finding a column 



density of (5.5±4.4) x 10^*^ cm'^, and Knez et al. (20091 



speculate that gas-phase OCS should be observable to- 
ward IRS 9 as indicated in some chemical models. We 
find an upper limit of iV(OCS) < 2.9 x 10^^ cm^^ in 
the gas phase, which is compatible with the Gibb et al. 
result to within errors; we therefore conclude t hat any 
PCS n ear IRS 9 is firmly sequestered as an ice. Tielens 
(20051 cites an OCS ice abundance of 0.1% relative to 
H2O ice toward IRS 9 given a water ice column density 
of 1 X 10^^ cm~^. That figure may be compared to a CO 
ice abundance of 10%. In interstellar ices, OCS is com- 
monly found in association with hydrogenated species 



such as H2O a nd CH3OH ( jPalumbo et al.j 119971 I gaF 
tois et al. 1 1999 ). Laboratory results reported oy Collings 
et al. 1 2004 ) indicate that OCS in an H2O mixture des 
orbs at ~150 K, implying that significant quantities of 
this molecule may be frozen out in cold envelopes. How- 
ever, C2H2 is seen toward both IRS 1 and IRS 9, for ex- 
ample, and according to the same study desorbs at only 
a slightly lower temperature. The relative columns of 
OCS and C2H2 may then reflect real differences in their 
abundances rather than their desorption temperatures. 
Alternately, OCS might only be abundant in the colder, 
outer envelope o f IRS 9, explaining its apparent absence 



in the gas phase. Knez et al. ( 2009 ) did not observe OCS 



in the gas phase toward IRS 1 and Gibb et aL] ([2C)04 ) only 



cite an upper limit of 1 x 10^^ cm"'^ for the column den- 
sity in the solid phase. Other sulfur-bearing molecules 
desorb from grains at a relatively low temperature of 
^70 K, well below the temperatures in both objects. It 
is difficult, then, to understan d why high- J lines of CS 
were observed to ward IRS 1 by Knez et a l. (2009) when 
recent models by Wakelam et alj ( l!01l[ ) underpredict the 
abundance of CS relative to other S- bearing molecules 
at temperatures above 100 K. The IRS 1 CS gas-phase 
column density is comparable to the solid-phase column 
density of OCS toward IRS 9, so if IRS 1 began with a 
comparable OCS abundance it may indicate a reaction 
route pre ference for forming CS a s opposed to, e.g., SO2. 



However, Wakelam et al. (2004) do not associate OCS 
and CS but rather assunie CS exists in the gas of the 
pre-evaporative phase. We did not search for CS or SO2 
toward IRS 9 on the basis of their presumed low abun- 
dances, and thus are unable to comment on the likelihood 
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of any OCS/CS association. 

The occurrence of isocyanic acid (HNCO) may be tied 
to the presence of the cyanate ion (OCN~), where OCN~ 
is formed on grains by irradiatio n of sohd HNCO in ices 
by U V in the presence of NH3 (van Broekhuizen et ah 
20041. Once evaporated from grains, OCiN~ may then 
undergo hydrogenation to become HNCO again. We 
speculate, however, that the hydrogenation of OCN~ 
takes place before desorption, since the charge an ion 
induces while on a grain w ould result in a relat ively high 
desorption temperature. Knez et al. (2009j) observed 



dozens of HNCO lines toward IKS 1 and mter a column 
density of 5.4 x 10^^ cm~^; this and the value of other 
gas parameters from their spectra led them to conclude 
that gas-phase HNCO may originate from evaporation 
of grain mantles by way of OCN~. A low column den- 
sity compared to the value expected from the observed 
abundance of OCN~ is interpreted as evidence that only 
a fraction of the available OCN~ desorbed from grains 
becomes HNCO, the remainder forming other molecules 
in gas-phase chemistry. Interpretation of observations is 
difficult, however, as Knez et al. point out that most 
chemical models do not currently include HNCO in their 
networks. All of the HNCO hues in the various TEXES 
settings we observed were non-detections. Our fitting 
program preferred column densities of (0.1±0.5) x 10^^ 
cm~^ and (0.9±1.6) x 10^^ cm~^ in the two presumed 
absorbing units, although a clear lack of detected lines 
make these values limits at best. We cannot, there- 
fore, confidently distinguish between scenarios in which 
HNCO has a very low abundance in all phases and those 
in which it is effectively locked up in ices and unavailable 
to spectroscopy at our wavelengths. 

5.3. Radiative Excitation and Radiative Transfer 

In order to explain the observed pattern of NH3 ab- 
sorption and emission in our spectra of IRS 9, we in- 
voke a radiative transfer effect that emerged from our 
modeling efforts. We begin with gas in a shell of thick- 
ness t around the continuum source with an NH3 density 
nNHs- The rotational levels {J,K) of the ground vibra- 
tional state are assumed to be populated thermally at 
temperature T. Electric dipole transitions to and from 
an upper vibrational level proceed according to the se- 
lection rule that AK = 0, AJ = ± 1 or 0. In ammonia 
there is also a doubling of each vibrational level into sym- 
metric and antisymmetric (s,a) inversion states, and the 
requirement that s ^ a. Transitions are labeled by the 
inversion symmetry of the lower level. This gives the 
permitted transitions a/s R{Ju-l,K), a/s Q{Ju,K) and 
a/s P{Ju+l,K). We obtain absorption coefficients a for 
each of these transitions from GEISA. 

The absorption optical depth corresponding to a is 



t{w) = a7V(NH3)0. 



(13) 



in which A'^(NH3) is the radial column density through 
the shell and 4)^ is the line shape function, taken to be a 
Gaussian. Molecules are excited to the upper state along 
all branches at a rate 



''up — 

P,Q,R 



Lw 
hew 



T{w)dw 



(14) 



is the incident continuum luminosity. The rate at which 
photons are emitted on a given branch is determined by 
the rate of transitions populating the upper level and the 
branching ratio of the Einstein A coefficients, e.g., 



Ap + A 



Q 



A, 



(15) 



for the P branch. We calculate values for the A coeffi- 
cients according to 



A^i = ^«(r)^e^'/'=^ 
y 9u 



(16) 



where A is the wavelength of the transition, Qr is the 
rotational partition function, gu is the statistical weight 
of the upper state, and Ei is the energy of the lower state. 

We consider the aP{A,K), aQ{3,K), and aR{2,K) lines 
that go to the J = 3 rotational level of the 51^2 vibra- 
tional level. We use the branching ratios to compute the 
expected emission flux in a given line I with wavenumber 
center w from 



F, 



Ap + Aq+ Ar 47rd2 z-.^ hew 



(17) 



To calculate the net flux in a line we also need the amount 
of absorption. We require the relative intensity of the 10 
fim dust emission feature at the frequency of each of the 
branches in order to compute the absorption line flux. 
We use mid-infrared spec tra of the you ng stars DI Cep 
and DK Tau obtained by Hanner et al. (1998) to which 
they fit models of optically-thin silicate dust emission. 
From these spectra we calculate relative intensities of 
0.52, 0.78 and 1.00 for the continuum at the 11.7, 10.7, 
and 9.7/im wavelengths of the aP, aQ, and aR lines for 
which J„ = 3, respectively. The intensities are normal- 
ized to the value on the aR branch as it is nearest to the 
peak of the 10 /zm silicate feature. The observed contin- 
uum flux near the P-branch lines, ~60 Jy, is turned into 
a luminosity at the source by multiplying by And^ , where 
d is the distance to NGC 7538 (2.65 kpc). We corrected 
for exti nction using the 10 ^ m optical depth measure- 



ment of Willner et al. 
2.2 magnitudes and A. 
9 



(|1982|), which implies Anjf^m ^ 
4.4 magnitudes toward IRS 



The extinction-corrected luminosity ■ 



per second assuming r <C 1 through the shell, where 



^1.5x 10^^ erg s^^ 

cm~^ at the position of the P-branch lines is scaled to 
the other branches according the values above. We as- 
sume throughout that the lines are sufficiently broad that 
they remain optically thin, leading to uncertainty in the 
covering factor. Only the product of column density and 
covering factor is meaningful under such circumstances, 
so we used the sum of this product for the two NH3 com- 
ponents rather than the sum of the NH3 column densities 
themselves in calculating P; . We considered the possibil- 
ity that the optically thin case might result in uncertain- 
ties propagating to the line fluxes, and to this end we 
carried out additional runs of FITSPEC in which we de- 
liberately held the value of the covering factor fixed equal 
to 1. This resulted in a variation of the product of column 
density and covering factor of < 20% compared to allow- 
ing the covering factor to vary. We therefore retained 
the assumption of the optically thin case, and opted to 
use the lower value for the product of column density 
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TABLE 6 

Predicted fluxes for the aP, aQ, aR 

LINES OF U2 TRANSITIONS OF NH3 (Ju 
= 3) AND MEASURED FLUXES FOR THE 
CORRESPONDING aP-BRANCH LINES 
FROM TEXES SPECTRA. 







Predicted 




Observed 


K 






FaR 









5.55 




-1.94 


1.80 


1 


2.64 


0.03 


-0.95 


1.35 


2 


2.03 


0.02 


-0.73 




3 


2.11 


-0.87 




0.90 



Note. — All fluxes are given in units 
of 10^** erg cm~^. Positive flux val- 
ues indicate net emission while negative 
values indicate net absorption in a given 
line. 



and covering factor of 8.52 x 10 cm in place of N in 
Equation 
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For a{T), we chose an excitation temper- 
ature of 2'DD' K, consistent with the FITSPEC results for 
the -61.6 km s^^ component. Finahy, the assumed ex- 
tinction was re-introduced at the end of the calculation 
to realistically predict the line fluxes observed at Earth. 

The line fluxes we compute are given in Table [6] The 
predicted fluxes are the differences of the strengths of the 
predicted emission and absorption components of each 
line and show an overall pattern of emission on the P 
branch, weak absorption or no lines on the Q branch, 
and absorption on the R branch. The P-branch fluxes 
compare favorably with the integrated line fluxes from 
our data implying that our explanation is at least plau- 
sible. This suggests that the dust and NH3 are in close 
proximity to one another, if not actually mixed. 

In the case of acetylene, we observed lines of the Q 
and R branches at 12.5-13.7 fim and the i^i + iy^ P branch 
near 7.5 /um. The 1/5 Q- branch and 1^4 + 1^5 P-branch lines 
are of low S/N, due to telluric interference, but clearly 
lead to higher C2H2 column densities than the 1^5 R- 
branch lines when the observations are interpreted with 
the pure absorption model. The 1/4 + 1/5 lines appear to 
require at least an order of magnitude more C2H2 than 
do the 1^5 R- branch lines. The radiative transfer effects 
discussed for NH3 may also be relevant for C2H2, and 
may explain a part of the discrepancy between the C2H2 
bands. But the bands are farther from the silicate feature 
and the wavelength separation between the branches of 
the C2H2 bands is smaller than for NH3, so the effects 
should be smaller. The greater column density derived 
from the + 1^5 lines could instead be due to the smaller 
extinction at 7.5 fim allowing a greater column of gas 
to be probe d. This explanation was suggested by |Evansj 
et al. ( 1991 1 to explain a similar effect seen toward Orion 
lKc2, but it is difficult to understand how it could explain 
more than a factor ~ 2 discrepancy, given the ratio of 
extinction at the different wavelengths. 

We propose a different explanation for the variation be- 
tween the column densities derived from the C2H2 and 
1^4 + 1^5 lines. Vibrational excitation caused by a band 
photon is almost always followed by emission of another 
1^5 band photon, although possibly in a different branch. 
(Due to the C2H2 selection rules, P- or P-branch absorp- 
tion can be followed by either P- or P-branch emission. 



but Q-branch absorption must be followed by emission 
in the same branch.) However, 7.5 /ini 1^4 + 1^5 absorp- 
tion is most often followed by a 1^4 -f 1^5 — i>4 transition 
at ^ 13.7 fim to the 1^4 level. Since radiative decay from 
the 1^4 level is forbidden, the 1^4 excitation energy is then 
eventually lost in a collision. The result of this series of 
transitions is that there is little 1^4 + emission cancel- 
ing the V4 + absorption, and these lines are probably 
reasonably well treated with the pure absorption model 
that we used. In contrast, the absorption lines may 
be largely cancelled by emission, leading to an under- 
estimate of the C2H2 column density. The amount of 
re-emission into our beam is very dependent on the (un- 
known) source geometry, which determines whether our 
line of sight has more or less than the average for the 
gas surrounding the source. Assuming the V4 -\- ^5 lines 
can be treated with a pure absorption model to give a 
fair measure of the abundance along our line of sight, the 
required correction to the column densities derived from 
the lines is approximately a factor of 10. The required 
correction to the column densities of our other observed 
molecules may be similar. 

5.4. Abundance Trends and Relative Ages 



In Figure 13 we present a graphical comparison of the 
column densities of C2H2, HCN, NH3, HNCO, and CH4 
we obtain tow ard IR S 9 with those derived for IRS 1 
( [2009t . 



Knez et al. 



Column densities for each of the 
two components with similar LSR velocities (-57 km s~^ 
and -61 km s~^) toward each object are shown. We 
reach rather different conclusions for the two velocity 
components; IRS 1 appears to be clearly enhanced in all 
molecules except for HNCO at -57 km s~^, whereas the 
column densities arc essentially the same to within errors 
at -61 km s^^. HNCO is the exception at the higher neg- 
ative velocity, but we again stress that the cited values 
for IRS 9 are those preferred by our fitting program de- 
spite HNCO being an evident non-detection in our data. 
It is therefore not strictly reliable in comparison with 
the firm detection of HNCO toward IRS 1. Also, we 
do not have any evidence a priori that these two velocity 
components are kinematically related. Knez et al. (20091 
claim the two components seen toward IRS 1 are real but 
may not be separate; they may instead represent a non- 
Gaussian velocity distribution. If the latter situation is 
true, then such a velocity distribution may be systemic 
in other parts of NGC 7538. Regardless, the pattern that 
emerges from Figure [13] is suggestive of real differences 
in the gas-phase abundance of these molecules, which 
we attribute to differing thermal conditions in the envi- 
ronment of each object. If the molecular abundances in 
all phases were the same to begin with, the gas-phase 
abundances we see at the present time are indicative of 
a difference in the relative ages of IRS 9 and IRS 1. 

Several other lines of evidence also argue in favor of 
a younger age for IRS 9. First, a number of molecules 
(e.g., H2O, CH3OH, HCOOH, and OCN") are seen in 
the solid phase in the spectrum of IRS 9 that are either 
in the gas phase or apparently absent toward IRS 1. CO 
and CH4 are seen both in the solid and gas phases, to- 
ward IRS 9 although the gas-phase lines probably form 
in the shocked outflow material. Second, a fairly strong 
ionizing radiation field exists near IRS 1, whereas our 
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Fig. 13. — Log column densities for the main molecular species detected in this stu dy for the two prin cipal velocity components toward 
each target. IRS 9 values (this work) are shown as filled circles and IRS 1 values from |Knez et al!] ( |2009^ as filled squares. The symbols for 
a given molecule are slightly offset for clarity. 



non-det ectio n of the 12.8 fim line of [Ne II] discussed in 
Section |4.4| argues for a much weaker field toward IRS 
9. In addition, IRS 1 has an ionized component to its 
outflow, probably driven by a wind; IRS 9 shows no sim- 
ilar outflow component. Third, the chemistry of IRS 
1 resembles the "hot core" phase in which ices desorb 
from dust grains, leading to a short-lived (~ 10^-10^ yr) 
high-temperature gas-phase chemistr y and the synthesi s 
of complex organic molecular species ( van der Tak|2005| . 
Some species, such as HNCO, are clear detections toward 
IRS 1 but appear to be absent toward IRS 9; we cannot 
rule out an intrinsic abundance effect in explaining this 
difference, but a simpler explanation is that the IRS 9 
c ore is not vet sufficiently warm to evaporate it. 

Knez et al.| (|2009) claim an age of ~10^ yr for IRS 1 
based on their chemical models. However, there is sig- 
nificant uncertainty in that value. Abundances of the 
molecules they detected indicate an age closer to ~10'' 
yr, but for other molecules the chemical model predic- 
tions disagree with the observed values, suggesting ages 
in the range of 2 x lO'^ to 2 x 10® yr. They explain this 
large range as a result of insuffic iently detailed models. 

(2002) predict enhancements 
, HON, and CH4 for temper- 
atures above about 800 K at late times (> 10^ yr); the 
enhancements are similar to those observed to ward IRS 
1, but the temperatures derived by Knez et al.' (2009) 



The models of Doty et al 
in the abundances ot C2H2 



are 



considerably lower than this. The fractional abundance 
of HCN in the Doty et al. model for 200 K quickly turns 
down at times > 10^ yr whereas those of C2H2 and CH4 
increase through at least 10^-^ yr; A^(HCN) is in fact It is tempting to assume 



slightly higher toward IRS 1 compared to IRS 9 by a fac- 
tor of about 1.4. This may both support the argument 
that IRS 1 is the older object of the two and suggest an 
upper limit for its age, since the HCN abundance has not 
begun the predicted decline, but this conclusion is tenta- 
tive since our value of iV(HCN) for IRS 9 is determined 
by the singular detection of the 1^2 i?(16) line. Lastly, 
an additional age constraint from chemistry is provided 
by the non-detecti on of the molecule CCS by Hoffman fc| 
Seojin Kim (2011) in EVLA data, suggesting an age > 
lO'' yr. Given that we find IRS 9 to be in an evidently 
more primitive chemical state than IRS 1, the upper limit 
to its age is probably ~ 10^ yr. Detailed, robust chemi- 
cal modeling, using our observations as inputs, may help 
more firmly constrain this value. 

Our results, when considered in the context of the ear- 
lier work on IRS 1, qualitatively support a picture of an 
evolutionary sequ ence of cores in NGC 753 8 co nsistent 
with the models of Elmegreen fc Lada ( 1977[) and]Camp-j 
bell & Thompson ( 1984[ ) in which a wave of shock/ ion- 
ization fronts from a previous generation of stars caused 
the sequential collapse of the cores that became IRS 
1 and IRS 9, respectively. As the wave continues to 
propagate toward the southeast, it should initiate the 
collapse of cores found in progressively more primitive 
states. The 'IR-quiet' protostellar object NGC 7538 S is 
likely the next younge st object in the region after IRS 9 
|Pestalozzi et al.|2006[ ). Such a sequence, including NGC 
7538 IRS 11 among the y oungest objects, is supported b: ' 
the maser observations of 



Hutawarakorn & Cohen (2003). 
hat the wave has left the IRS 1 
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region already, but several newly-discovered s ubmillime- 
ter sou rces within 0.35 pc of IRS 1 reported by |Qiu et al. 
(2011), presumably harboring forming intermediate- or 
high-mass stars, may be younger than IRS 9 and contem- 
poraneous with NGC 7538 S. However, they may also be 
fragments, along with IRS 1, of a common progenitor. 
The star formation history of NGC 7538 may therefore 
be more complex than can be represented with a simple, 
monotonically-increasing series of protostellar ages. 

5.5. Spatial Resolution of IRS 9 

We checked whether IRS 9 is spatially resolved in 
TEXES data by obtaining spatial-spectral maps of IRS 
9 and n Cep, a mid-infrared- bright standard star that 
was used as a PSF reference for spatial deconvolution. To 
generate the maps, the TEXES slit was stepped across 
each object and spectral information gathered at each 
step, building up data cubes in which there are two spa- 
tial dimensions and one spectral dimension. Additional 
steps were included at both the beginning and end of the 
scans to sample the sky for removal. The data are han- 
dled by our reduction pipeline in a manner substantially 
similar to the method of processing the nodded observa- 
tions, treating a scan as a series of 2-D spectrograms. 

IRS 9 and /i Cep were observed with identical scan pa- 
rameters including scan direction, length, and step size, 
and were observed at comparable zenith angles. After 
reducing the raw data, different spectral windows were 
selected from which to form reconstructed images of each 
object by summing the data in the spectral direction. 
For IRS 9, we chose a window centered at about 745.88 
cm~^, near the position of the C2H2 1^5 R{G) line in the 
rest frame of the object. The window was 0.35 cm~^ 
wide. We selected this frequency range to search for spa- 
tial extent in the IRS 9 map as we noted absorption in 
the R{6) line in nod-mode spectra. For the deconvolution 
reference, we chose a nearby spectral region apparently 
free from telluric and stellar photospheric absorption cen- 
tered at 744.43 cm"^ and 0.25 cm"^ in width. The ob- 
served FWHM for each object was 0'.'67 for IRS 9 and 
0'.'63 for fj, Cep at each object's respective spectral set- 
ting. The deconvoluti on was performed using the maxi- 
mum entropy met hod ( [Bryan &: Skilling 1980 Nara yan fc 
Nityananda| 1986 1, with results shown in l^'igure L'he 
deconvolvecf image of IRS 9 in the C2II2 feature appears 
pointlike with some possible weak extension toward the 
NE. We do not believe this indicates that the envelope 
or other structure of IRS 9 is spatially resolved, and we 
carried out our other TEXES observations exclusively in 
nod mode. 

This approach to imaging and deconvolution should 
be interpreted cautiously. Poor seeing during a scan, for 
example, might distort the PSF on time scales shorter 
than the time required to scan across an object. How- 
ever, we felt sufficiently convinced that IRS 9 presents 
itself as a point source in the mid-infrared that we did 
no further scans. Given the distance to NGC 7538 and 
the spatial resolution of our observations, we conclude 
that spatial structure of IRS 9 is unresolved on scales < 

2000 AU. This limit is consistent with the inferred sizes 

i3 1 n5 



(Patel et al. 2005| Schreyer et al. 2006 


Rodn'guez et al. 


2007 


l-'ranco-Hernandez et al.|2009|lGa 


van-Madrid et al. 


2010 


1 and places a constraint on physical models of mas- 
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Fig. 14.— Contour plot of the TEXES scan map of NGC 7538 
IRS 9 at 745.8 cm~^ after maximum entropy deconvolution. The 
PSF reference for the deconvolution was a scan map of the mid- 
infrared standard fj, Cep made with the same parameters as the 
IRS 9 scan. North is up and east at left and contours are plotted 
at flux densities of 25, 50, 100, 150, 200, 250 and 300 Jy. The plot 
origin is approximately the location of the source peak intensity. 

sive protostars in NGC 7538, which often include disks 
thought to be hundreds of AU in size (e.g., Sandell et al. 



2009 



fiht to be hundreds c 
Surds et aL]|2011 |. 



5.6. A Structural Model Of IRS 9 

Based on the conclusions we reach from kinematic in- 
dications in the data, we compose a toy model of the 
IRS 9 system to show its essential structure and a guess 
at the viewing geometry. This model is shown in Fig- 
In this view, the central, perhaps accreting ob- 
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ject is shown in cross-section with a surrounding disk. 
The known bipolar outflow is indicated by the evacuated 
regions above and below the disk and by arrows showing 
the direction of the flow. A lighter shading above and 
below the disk represents the remaining envelope that 
has not been accreted or become part of the outflow. In 
order to account for seeing certain molecules under pre- 
sumed LTE conditions while not observing the exciting 
source directly, we suggest a "dust photosphere" at a ra- 
dius from the source at which the emission from the dust 
becomes optically thick, as indicated by dashed lines in 
the figure. It is this effective source of radiation against 
which we see most of our absorption lines. Since we do 
not resolve spatial structure in the molecular absorption 
region above the dust photosphere, we interpret the size 
scale constraint discussed in the previous section (2000 
AU) as an upper limit to the radius of the dust photo- 
sphere itself. This value also represents an upper limit of 
the size of the disk, because a larger, optically-thick disk 
should have been detectable by virtue of its spatial ex- 
tent. Our data do not permit a constraint to be placed 
on the size of the surrounding envelope or the outfiow 
cavities. 

The inclusion of a disk in this picture is tentative, given 
that the lifetime of disks around massive stars is thought 
to be particularly short, of order 10~^ the lifetime of the 
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< 2000 AU 



Fig. 15. — A cartoon rendering for a proposed model of NGC 
7538 IRS 9 based on TEXES data. A central, embedded source 
(or sources) has at least partially cleared cavities in a bipolar out- 
flow while significant dust remains in a spherical distribution (light 
gray shading) . An inner disk may exist (dark gray shading) . To ac- 
count for our observations, the line of sight samples both scattered 
light in the outflow as well as a "dust photosphere" (dashed lines) 
where the dust becomes opaque at a mid-infrared wavelengths. The 
inferred viewing angle is shown, and a scale bar is provided based 
on the constraints discussed in Section [5.5[ Note that the envelope 
and outflow cavities as depicted are not shown to this scale. 



host star (Mann & Williams 2009). The evolutionary 
state of iRS 9 may be sutticiently early in its history, 
well be fore the embedded protostar is visible, to retain 
a disk. Knez et al. (20091 argue for a disk around IRS 
1 whose surtace is observed at near-grazing incidence; 
a sightline through the disk atmosphere is proposed to 
explain their observations. 

We propose a viewing angle based on the following. 
First, the broad CO lines with their large, Doppler- 
shifted velocity separations suggest formation in the 
bipolar outflow. Second, the apparent radiative "pump- 
ing" of the NH3 transitions, discussed in Section |5.3[ 
requires continuum light from the direction of the dust 
photosphere. Third, the other molecular species we ob- 
served in absorption also require a background contin- 
uum source, though the range of excitation conditions 
we derive is insufhcient to definitively locate them in the 
envelope or the turbulent edge of the outflow cavity. The 
C2H2 combination mode may be an exception and could 
shock-excited in the outflow. Consequently, we place the 
viewing angle such that the observer is looking nearly 
along the edge of the outflow cavity and sees light from 
both sources. This model appears to account for all as- 
pects of our observations in at least a qualitative way. 
It is a lso consistent with the conclusion of Sandci ret al.l 
(2005), based on HCO'^ J = 2— >1 observations show- 
ing deep, redshifted self-absorption, that the molecular 
outflow is viewed nearly pole-on. 

6. SUMMARY 



We have presented the method and results of a study of 
the embedded high-mass protostellar object NGC 7538 
IRS 9 and compared our findings to a similar investiga- 
tion of the related object NGC 7538 IRS 1. We obtained 
high resolution, mid-infrared spectra of these objects in 
46 ro-vibrational transitions of the fundamental bands 
of the molecules C2H2, CH4, HCN, NH3 and CO and 
a number of their isotopologues. We also detected two 
lines of the -\- 1^5 combination mode of C2H2. From 
these observations we draw some broad conclusions. 

1. IRS 9 appears to be a spatially- unresolved object 
on a scale of -2000 AU. 

2. With the exception of CO, whose lines are satu- 
rated, we did not observe ^'^C isotopologues of any or- 
ganic molecules. This may indicate a particularly large 
value of the ^^C/^^G ratio, or the ^'^C isotopologue lines 
were simply below our threshold of detection. 

3. There is no discernible, consistent trend in the abun- 
dance variations of C2II2, HCN, NII3, and CH4 with 
respect to either CO or II2 between IRS 9 and IRS 1. 
Column density variations between the two objects are 
an order of magnitude or less in each case. 

4. The observation of gas-phase CS absorption toward 
IRS 1 and the evident non-detection of sulfur-bearing 
species like OCS in the gas phase toward IRS 9 may 
reflect a real S abundance variation between the two ob- 
jects but can also be explained by different thermal con- 
ditions in each case. A similar abundance discrepancy 
involving HNCO may indicate sufficiently low tempera- 
tures in the IRS 9 envelope to retain it as an ice and 
explai n its detection in the solid phase by |Gibb et al.| 
( 2004 1 . If the abundances in all phases are actually com- 
parable between the two sources, a relatively weak heat- 
ing source is implied in the case of IRS 9. While this is at 
odds with the somewhat higher inferred abundances of 
the other molecules toward the warmer environment of 
IRS 1, the suggested sightline in that direction through a 
disk atmosphere subjects those molecules to direct irradi- 
ation by protostellar UV. This may drive photochemistry 
that lowers abundances of certain molecules as they are 
converted to other kinds. 

5. However, the relative abundances of various molec- 
ular species with respect to CO and H2 can also be ex- 
plained by variations in the CO/H2 ratio or the gas-to- 
dust ratio between the two objects. Choosing between 
these interpretations would require better determining 
the CO and H2 abundances toward both objects. 

6. Lines of many species observed toward IRS 9 have 
higher Doppler b values than those toward IRS 1, result- 
ing in lines toward IRS 9 being less saturated than those 
toward IRS 1. 

7. Given the non-detection of Ne II in our data and the 
continuum flux density at the expected wavenumber of 
the line (70±10 Jy), we find a corresponding upper limit 
of —10^^ photons s~^ for the Lyman continuum flux in 
IRS 9. This value is well below both the expected flux of 
a ZAMS star with the same luminosity and the ionizing 
flux of IRS 1 (> lO'^s photons s^^). 

8. The observed pattern of emission and absorption 
components of ammonia lines in our spectra can be ex- 
plained by a non-LTE radiative transfer effect involving 
"pumping" of lines on the i?-branch by light from the 9.7 
/im silicate dust feature. Predictions of the strength of 
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this effect compare reasonably well with observed fluxes 
in P-branch lines. 

9. Our data are consistent with a simple model of the 
IRS 9 system in which the sightline probes both a high- 
speed outflow and a quiescent envelope illuminated from 
behind by a "dust photosphere" such that line formation 
depth is wavelength-dependent. The possibility of seeing 
the turbulent flow along the wall of the outflow cavity 
is supported by the broad, blueshifted wing of the C2H2 

+ Vb combination band that may arise in a shock. A 
disk is included in the model given the system's inferred 
young age and the presence of a kinematically-indicated 
disk toward the (presumably) older IRS 1. 

10. Observed differences in the abundance of various 
molecular species toward IRS 9 and IRS 1 imply an upper 
limit to the age of IRS 9 of ^ 10'' yr, whereas IRS 1 
may be an order of magnitude older. This conclusion 
is consistent with th e trigg ered star formation models of 
Elmegrecn & Lad aJ (|1977 j as apphed to NGC 7538 by 
Campbell & Thompson 1 1984 ) and others. 
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